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In Hofava and Einstein-aether theories of modified gravity, in spite of the violation of Lorentz 
invariance, spherically symmetric stationary black hole solutions possess an inner universal horizon 
which separates field configurations into two disconnected classes. We compute the late time radia¬ 
tion emitted by a dispersive field propagating in such backgrounds. We fix the initial conditions on 
stationary modes by considering a regular collapsing geometry, and imposing that the state inside 
the infalling shell is vacuum. We find that the mode pasting across the shell is adiabatic at late time 
(large inside frequencies). This implies that large black holes emit a thermal flux with a temperature 
fixed by the surface gravity of the Killing horizon. In turn, this suggests that the universal horizon 
should play no role in the thermodynamical properties of these black holes. 
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I. INTRODUCTION 

The laws of black hole thermodynamics are firmly established in Lorentz invariant theories, and they play a crucial 
role in our understanding of black hole physics [1] . In particular, the entropy and the temperature are governed by 
the area and surface gravity of the event horizon. In Lorentz violating theories, the status of these laws is unclear 
because essential aspects are no longer present [2-7]. For instance, the thermality of the Hawking flux is inevitably 
lost in the presence of high frequency dispersion, although it is approximatively recovered for large black holes, i.e., 
when the surface gravity k is much smaller than the UV scale A setting the high frequency dispersion [8]. 

The origin of the difficulties can be traced to the fact that the event horizon no longer separates the outgoing field 
configurations into two disconnected classes. In fact, when the dispersion is superluminal, it can be crossed by outgoing 
radiation. However, it was recently discovered that in some theories of modihed gravity such as Hofava gravity [9-11] 
and Einstein-aether [12-15], spherically symmetric black hole solutions possess a second inner horizon. This horizon, 
named universal, cannot be crossed by outgoing configurations, even for superluminal dispersion relations which allow 
for arbitrarily large group velocities. (The difficulty mentioned in [2] is thus evaded.) Following this discovery, it 
has been argued that the universal horizon should play a key role in the thermodynamics of such black holes. Even 
though they seem to obey a first law [16, 17], a key question concerns the temperature of the Hawking radiation they 
emit. Would it be essentially governed by the (higher) surface gravity of the universal horizon, or would it still be 
fixed by the surface gravity k of the Killing horizon? 

Two recent works concluded that the universal horizon emits a steady radiation with properties governed by its 
surface gravity. Because of the complicated nature of the field propagation near that horizon, this conclusion was 
indirectly obtained, in [18], by making use of a “tunneling method”, and, in [19], by analyzing the characteristics of 
the radiation field. In the present paper, we reexamine this question by performing a direct calculation and reach the 
opposite conclusion that no radiation is emitted from the universal horizon at late time. 

We proceed as follows. As in the original derivation of Hawking [20], we identify the boundary conditions on the 
outgoing modes in the near vicinity of the universal horizon by considering a simple collapsing shell geometry, and 
by assuming that the state of the field is vacuum inside. We then compute the mode mixing across the shell between 
inside modes propagating outwards and outside stationary modes with a fixed Killing frequency t/jx. The late time 
behavior is obtained by sending the inside frequency w —t oo. In this limit, we show that the scattering coefficients 
involving modes with opposite norms vanish. This result can be understood from the fact that the modes are 
accurately described by their WKB approximation in the immediate vicinity of the universal horizon. In other words, 
the pasting across the shell is adiabatic in the limit w —?► oo. Hence, for large outgoing radial momenta, the state of 
the field outside the shell is the usual vacuum, as explained in [21]. 

It then remains to propagate these high momentum dispersive modes from the universal horizon till spatial infinity. 
This propagation has already been studied in detail; see [22] for a recent update. It establishes that large black 
holes emit a stationary flux which is (nearly) thermal, and with a temperature approximatively given by the standard 
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relativistic value. In other words, the robustness of the Hawking process, i.e. its insensitivity to high frequency 
dispersion which was first established in [23], is now extended to black holes with a universal horizon. 

From this it is tempting to conclude that the laws of black hole thermodynamics should also be robust, and 
they should involve the properties of the Killing horizon. This conclusion is reinforced by the fact that the field 
configurations propagating on either side of a universal horizon come from two disconnected Cauchy surfaces, and are 
highly blueshifted. Hence, it seems that no Hadamard condition of regularity [7] could be satisfied on the universal 
horizon. This raises the question of the fate of the universal horizon; see [6]. This difficult question shall not be 
discussed in the present work. 

Appendix A gives the details of the calculation which is summarized in the main text. In Appendix B we compare 
our model with previously-studied dispersive ones without a universal horizon, and show the role of the acceleration 
of the preferred frame. Appendix C shows the results of numerical simulations confirming the approximately thermal 
character of the emission at infinity governed by the surface gravity of the Killing horizon. 


II. MASSLESS RELATIVISTIC SCALAR FIELD IN A COLLAPSING SHELL GEOMETRY 

In this section, we briefly review the computation of the Hawking radiation emitted at late time in a collapsing 
geometry [20] . Although these concepts are well known, we present them in a way which prepares the more involved 
calculation of the late time flux when dealing with a dispersive field in the presence of a universal horizon. As explained 
in the Introduction, we shall use a direct calculation which consists of pasting the modes across the infalling shell. 
We closely follow the derivation of [24]. 

For simplicity, we consider an infalling spherically symmetric lightlike thin shell. In this case, it is particularly 
appropriate to work with advanced Eddington-Finkelstein (EF) coordinates f,r, where v is the advanced null time. 
At fixed r, one has dv/dts = 1, where ts is the usual Schwarzschild time. Hence, outside the shell, the stationary 
Killing field is simply On both sides of the shell taken to be u = 4M, the line element reads 

ds^ = _ 2dvdr — -I- sm{9)'^d(p^) . (1) 

where M{v) = Q{v — AM)M. These coordinates cover the entire space-time, shown in the right panel of Fig. 1. On 
the left panel, the infalling and outgoing null radial geodesics are represented in the (v — r, r) plane. One clearly sees 
that the Killing horizon (where the norm vanishes) divides the outgoing geodesics into two separate classes. 

We work in Planck units: c = h = G = I. 

Let $ be a massless real scalar field with the action 

.5 = (2) 

We define ip = and consider radial solutions independent of {9, (p). Inside the shell, for v < AM, we introduce the 
null outgoing (affine) coordinate U = v — 2r. Outside the shell and for r > 2M, we introduce the null coordinate 

u = V — 2r)f G (— oo,oo), (3) 

where r’^ = r + 2Mln |r/2M — I| is the usual tortoise coordinate, which diverges on the Killing horizon. To cover 
the region inside the Killing horizon, one needs another coordinate ul = —{v — 2r^). ^ The field equation then reads 

dudyip = 0, V < AM, , , 

(5„a„ + (l-^)^)V- = 0, v>AM. 

For simplicity, we neglect the potential engendering the grey body factor and work with the conformally invariant 
equations dudyip = dydyip = 0. The solutions can be decomposed as 

ip{u,v) = ip'^{u)+ip'"{v), (5) 

and similarly for v < AM with u replaced by U. The infalling v sector and the outgoing u sector completely decouple. 
Moreover, the v modes ip'" are regular across the horizon and play no role in the Hawking effect. We thus consider 
only the u modes, and, to lighten the notations, we no longer write the upper index u on outgoing modes. 


^ This sign guarantees that dU/du^ is positive. As we shall see in Section III, a similar sign must be taken when studying a dispersive 
field on both sides of a universal horizon. 
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FIG. 1. (Left panel) Null radial geodesics in the {v — r,r) plane in units of M. The coordinate v — r coincides with the 
Minkowski time T inside the mass shell, and with the Schwarzschild time ts for r —>■ oo. The solid black lines are null radial 
geodesics which are reflected on r = 0. The dashed green line represents the trajectory of the null shell v = 4M. The blue, 
dot-dashed one shows the Killing horizon at r = 2M outside the shell. The dotted purple line shows the locus r = M, v > AM 
which will play a crucial role in Section III. The wavy line shows the singularity, located at r = 0, u > 4M. (Right panel) 
Penrose-Carter diagram of the collapsing shell geometry. The vertical line corresponds to r = 0, u < 4M. X~ corresponds to 
u,U ^ — 00 , and to u —>■ oo. 


To compute the global solutions, we need the matching conditions across the null shell. In the present case, ip is 
continuous along v = AM. Hence ipinside{U) = ipoutsideiu{U)), where the relation between null coordinates is 


u{U) = U-AM\n 



( 6 ) 


for r > 2M {U < 0). For r < 2M {U > 0), one has ul{U) = —u{\U\). 

To obtain the Hawking flux one needs to relate the in modes (p™ characterizing the vacuum inside the shell, to 
the out modes characterizing the asymptotic outgoing quanta with Killing frequency A. In the internal region, a 
complete orthonormal basis of positive-norm modes is provided by the plane waves 




e 


— iujU 


2^'KLd ’ 


( 7 ) 


where w € K’*' is the inside frequency id/y. In the external region, the (positive-norm) stationary modes for r > 2M, 
are 


= 0(r- 2M) 


^—i\u 


2^/PP\' 


A e 


( 8 ) 


A similar equation defines (P^^\ul) in the trapped region, for r < 2M. The modes and their complex 

conjugate form a complete orthonormal basis. One easily verifies that the conserved scalar product for the u modes 
can be written as 


du {iplduip 2 - 'ip 2 duipi) ■ 


(V'l,V'2) = * 


— oo 


( 9 ) 
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The Bogoliubov coefficients encoding the Hawking flux are then given by the overlaps between the two sets of 
modes: 


/3A,. = ((<^rr,c)- (10) 


Using u{U) of Eq. ( 6 ), they can be computed explicitly; see [24] for details. The late time behavior is obtained by 
sending the inside frequency uj —>■ oo. In this limit, one recovers the standard thermal result 


^ SttMA 

Q!A,w 


( 11 ) 


To prepare for the forthcoming analysis, it is instructive to compute the Bogoliubov coefficients by the saddle point 
method [25, 26]. For the ax,uj coefficient, when w » A, i.e., at late time, the location of the saddle is given by 


\ = U) 




( 12 ) 


where mq is a constant which drops out of the late time flux. (In the present model, uq vanishes.) From this equation 
we recover the time-dependent redshift relating w, the large frequency emitted from the collapsing star, to A, the 
frequency received at infinity and measured using the proper time of an observer at rest. In particular, we recover the 
characteristic exponential law governed by the surface gravity At = i9„ In C/ (u) = . Had we considered a collapsing 

shell following a (regular) infalling timelike curve, Eq. (12) would still have been obtained at late u — uq time. 

This is the kinematical root of the universality of Hawking radiation in relativistic theories. Indeed, when studying 
the /3 a, w coefficient, one finds that the saddle point is now located at A = . When taking into account 

the fact that the integration contour should be deformed in the lower it-complex plane, one finds that Us.p. has an 
imaginary part 9(its.p.) = —tt/ac, whereas its real part in unchanged. This gives a relative factor exp(— ttA/ac) with 
respect to the aA,^ coefficient. Upon squaring their ratio, we recover Eq. (11). We also recover here that the Hawking 
temperature k/2tt is fixed by the late time exponential decay rate entering Eq. (12). We finally notice that the 
stationarity of the flux is nontrivial. It follows from the fact that the ratio of Eq. (11) is independent of w, and from 
the fact that |/3 a,i.jP oc l/w for w —oo [26]. 


III. EMISSION FROM A UNIVERSAL HORIZON 
A. The model 


We aim to compute the late time radiation of a dispersive field propagating in a collapsing geometry. In principle, 
the radiation and the background fields should both obey the field equations of some extended theory of gravity, such 
as Hofava-Lifschitz gravity [9] or Einstein-aether theory [12, 13]. Since our aim is to study the radiation rather than 
the collapse, the latter shall be described by a simplified model. At the end of the calculations, we shall argue that 
our results do not qualitatively rely on the particular model we use. 

For reasons of simplicity, we assume that the collapsing object is a null thin shell, and that the external geometry 
is still Schwarzschild. In this case, the metric is again given by Eq. (1), and the Penrose diagram of Fig. 1 still covers 
the whole space-time. To describe the (unit time like) *ther field in the external region outside the shell, we 
adopt the solution of [16] (also used in [19]) with C 123 = 0, rg = 2M, and = 0. The Killing horizon is still at 
r = 2M, whereas the universal horizon, where = 0, is located at r = M. Inside the shell, we assume that the 

*ther field is at rest. To our knowledge, this configuration has not been shown to be a solution of the field equations. 
However, as explained in Subsection HIE, small deviations from this configuration should not significantly modify 
our conclusions. 

In EE coordinates, on both sides of the shell, the aether field u^, and its orthogonal spacelike unit field are given 

by 


u^dp. = dy 
s'^df, = dy 


MWg, 

r 



(13) 


where M{v) = Q{v — AM)M. We introduce the “preferred” coordinates t,X by imposing that oc dt and 

= sgn(r — M)dx- Their precise definition is given in Appendix A 1 . In these coordinates, the metric takes the 
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Painleve-Gullstrand form: 

where 


ds^ = c^dt^ - {dX - Vdt) 
M{v) 

M{v) 


V = -- 


c = \K^^u^\ = 


1 - 


At fixed t, outside the shell, V and c only depend on X. We notice that 


Ufj_dx^ = cdt. 


(14) 


(15) 

(16) 


The factor c ensures that dt is a total differential. Moreover, as explained in Appendix B, c is constant when the 
aether field is geodesic. Here we work with an accelerated aether, which is a necessary condition to have a universal 
horizon. Importantly, c vanishes on the universal horizon. ^ In fact, the novelties of the present situation with respect 
to the standard case studied in [8] only arise from the vanishing of c, and the associated divergence of the dispersive 
scale A/c. 

In Fig. 2 we show the lines of constant preferred time, and the direction of the aether field in the u,r plane. 
The coordinate t is discontinuous across the shell trajectory, as was the null coordinate u in the former section. As 
in the relativistic case, outside the shell we must use two coordinates t and now on either side of the universal 
horizon. The inside coordinate T evaluated along the shell, at u = AM~, is a monotonically increasing function of 
both t{v = 4M+,r) for r > M and of tL{v = 4M+,r) for r < M. So, the foliation of the entire space-time by the 
inside coordinate T is globally defined and monotonic. 

We consider a real massless dispersive field $ with a superluminal dispersion relation. Its action is given by Eq. (2) 
supplemented by a term quartic in derivatives: 


S = 



(a^$) (a'^d-) 


^ (V^ (Vp (/i^"V.$))l , 


(17) 


where Vp is the covariant derivative and — u^u’' is the projector on the hyperplane orthogonal to . The 

dispersive momentum scale is given by A. The field equation reads 

^ $ = 0. (18) 

Using a (1 -f l)-dimensional approximation, Eq. (18) reduces to 


[dt + dxV] ^ [dt + Vdx] - dxcdx + ^dxcdx^dxcdx 
c c 


V^ = o, 


(19) 


when working in the preferred coordinates. Since this (self-adjoint) equation is second order in dt, the Hamilton 
structure of the theory is fully preserved. In particular, the conserved scalar product has the standard form 


(^ilV’2) = i 


dX {^1^2 - nIV' 2 ) , 


( 20 ) 


where H = u^d^tj) = {dtip + Vdxip) /c is the momentum conjugated to '0. For more details; see Appendix A 2. 
The Hamilton-Jacobi equation associated with Eq. (19) is 


= {X- V{X)P f = c{X f 
We introduce the Killing frequency A, the preferred frequency H, and the preferred momentum P: 


pA 

A2 


A = -K^d^S = -dtS, 
n = -c(A) uf^dt^S = A - V{X)P, 
P = S^df^s = dxS. 


( 21 ) 


( 22 ) 

(23) 

(24) 


In these equations S should be conceived as the action of a point particle; see [21, 22, 29]. As explained in these 
works, S governs the WKB approximation of the solutions of Eq. (19). Notice that Eq. (22) only applies outside the 
shell, whereas all the other equations make sense on both sides. 


^ In an analogue gravity perspective [27, 28], to reproduce such a situation one needs a medium in which the group velocity of low- 
frequency waves vanishes at a point. From Eq. (21), we see that the effective dispersive scale A/c must be divergent at the point where 
c ^ 0. It would be interesting to find media which could approximatively reproduce this behavior. 
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FIG. 2. In this figure we show the lines of constant preferred time for the collapsing geometry in the plane {v — r)/M,r/M. 
The dashed line represents the trajectory of the null shell v = 4M, and arrows show the direction of the aether field u^. Notice 
that the external preferred time t diverges on the universal horizon r = M, v > 4M, whereas the internal time T, which is 
equal to v — r inside the shell, covers the entire space-time. 


B. The modes and their characteristics 

To compute the late time radiation one should identify the various solutions of Eq. (19), and understand their 
behavior. In the presence of dispersion, one loses the neat separation of null geodesics into the outgoing u ones, 
and the infalling v ones. In what follows, we call P“ (P") the roots of the dispersion relation which have a positive 
(negative) group velocity in the frame at rest with respect to the “fluid” of velocity V; see [8]. Similarly, the 
corresponding modes will also carry the upper index u ot v. 


1. The in and out asymptotic modes 


In the internal region v < 4M, the situation is particularly simple. Since the velocity field V vanishes, the preferred 
frequency is a; = —BtS, and the dispersion relation Eq. (21) becomes 


= p2 


p4 

A2- 


(25) 


This relation is shown in the left panel of Fig. 3. At fixed w, the positive frequency modes with wave vectors P“(a;) > 0 
and P'"{uj) < 0 define the two in modes and They both have a positive norm, which can easily be set to 

unity through a normalization factor. The mode is the dispersive version of the relativistic in-mode of Eq. (7). 
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FIG. 3. (Left panel) Dispersion relation in the internal region, where the preferred frame is at rest, in the D,P plane. The 
solid line shows lj = Q versus P for the positive-norm modes. The dashed line corresponds to negative D, i.e., negative-norm 
modes. The intersections with a line of fixed a; > 0 (dotted line) give the two solutions and P”. Right: Dispersion relation 
in the “superluminal” region for M < r < 2M in the A, P plane. The two additional roots on the negative D u branch are 
clearly visible. 


Outside the shell, for v > AM, at fixed Killing frequency A > 0, the situation is more complicated as the number 
of real roots depends on r. Outside the Killing horizon, for r > 2M, one has c > \V\. So, Eq. (21) possesses two real 
roots P“(A) > 0 and P"(A) < 0, which describe outgoing and infalling particles, respectively. The WKB expression 
for the corresponding stationary modes [the solutions of Eq. (19)] is 

exp(-i(At- p(^^A,X')dX')) 

^W (t,X)p, - \ ^ , (26) 

47r^|fl(A,PW)/(c(X)5vPW)| 

where P^’''>{X,X) is a real solution of Eq. (21) at a fixed A, and f2(A, P^®)) the corresponding preferred frequency. 
These WKB modes generalize the expressions of [22, 30] in that c is no longer a constant. Using Eq. (20), one 
easily verifies that they have a unit norm. One also verifies that the group velocity along the «th characteristic is 
dX^'^'>/dt = 1/i9aP^*^- When considered far away from the black hole, r/(2M) 3> 1, the u-WKB mode is the dispersive 
version of the relativistic out-mode of Eq. (8). 

From this analysis, we see that there is no ambiguity to dehne the asymptotic behavior of the in and out modes, 
solutions of Eq. (19). As before, these two sets encode the black hole radiation through the overlaps of Eq. (10). To 
be able to compute these overlaps, we need to construct the globally defined modes. To this end, we must study both 

(i) 

the behavior of ^ near the horizon and the third kind of stationary modes which propagate in this region. 


2. Near horizon modes 

Inside the Killing horizon but outside the universal horizon, for M < r < 2M, one has c < jUj. As can be seen 
from the right panel of Fig. 3, one recovers the two roots P“(A) > 0 and P®’(A) < 0 we just described. One notices 
that the u root P“(A) has been significantly blueshifted, whereas the infalling root P®'(A) hardly changed. Locally, in 
the WKB approximation, the corresponding modes are again given by Eq. (26). 

In addition, below a certain critical frequency Ac that depends on c and V, we have two new real roots we call 
where the arrow indicates the sign of the group velocity given by l/d\P. (The minus signs in 
front of these roots and A come from the fact that they have a negative preferred frequency 12 for A > 0. Hence, for 
A = — [A], the mirror image roots, Pi*|)[[^^ and have a positive 12.) Since 12 < 0, the WKB modes associated 

with these roots have a negative norm [22]. We call the right-moving one and the left-moving one , 

so that the modes without complex conjugation have a positive norm. Both of them carry a negative Killing energy 
—A. Using Eqs. Eq. (26) and Eq. (20), one easily verifies that and have a negative unit norm 

within the WKB approximation. As we shall see they describe the negative energy partners trapped inside the Killing 
horizon before and after their turning point, respectively. To summarize the situation, it is appropriate to represent 
the characteristics of the three types of modes. We proceed as in [21, 22]. 
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FIG. 4. Characteristics in a Schwarzschild stationary geometry for A = 10“^A (left) and A = A (right). The arrows indicate 
the direction of increasing preferred time along each characteristic. Solid lines correspond to positive-norm modes, and dashed 
ones to negative-norm modes. For r > M, each characteristic is named by the corresponding mode. The green dashed line 
corresponds to an extra v mode confined in r < M, as discussed in footnote 3. In this footnote, we also explain that the 
infalling v mode (described by the bine line) possesses a turning point inside the universal horizon. The mode corresponding 
to the orange, dashed line is the high momentum WKB mode before the turning point, and the low momentum mode 

(V’-’a”)* after it. 


3. The characteristics 

As said above, the characteristics are solutions of the Hamilton-Jacoby equation Since the frequency 

is a constant of motion on each side of the mass shell, they can then be computed straightforwardly. In Fig. 4, they 
are shown in the external region v > AM for a small value of |A|/A = 0.01 (left panel) and a moderate one |A|/A = 1 
(right panel). The solid lines correspond to positive energy solutions while the dashed ones correspond to negative 
energy solutions. 

The infalling n-like characteristics corresponding to V’a (in blue) approach the universal horizon from infinity and 
cross it at a finite value of v. (When sending A —>■ 0 they asymptote to null infalling geodesics v = cst.) As their wave 
vectors are finite for r M+, these characteristics will play no role in the sequel. As in the relativistic case, the v 
modes act as spectators in the Hawking effect. ^ 

The M-like characteristics with positive energy (in red), corresponding to the WKB modes emerge from the 
universal horizon from its right (r > M) at early times. When t increases, the momentum is redshifted while r 
increases. At a finite time, the characteristics cross the Killing horizon, and go to infinity as t > oo (almost along 
null outgoing geodesics when A/A ^ 1). 

The third characteristics (orange, dashed line) describe the trajectories followed by the negative-energy partners. 
For t —>■ —oo, they also emerge from r = M~^. However, when increasing t they have a turning point inside the Killing 
horizon, after which they move towards the universal horizon, smoothly cross it, and hit the singularity at r = 0 at 
finite values of v and Before the turning point, they are described by the WKB mode [tp-’^) ; after the 
turning point by (l//.’/”)*. 

It is important to notice that the only novel aspect with respect to the standard dispersive case (treated in full 
detail in [22]) concerns the behavior near the universal horizon. To clarify these new aspects, we represent in Fig. 5 
the global structure of the characteristics in the collapsing mass shell geometry. 


® Interestingly, r)-like characteristics have a turning point inside the universal horizon r < M. (The presence of the turning point may be 
understood from the fact that, close to r = 0, \V\ and c go to infinity but |y|/c goes to 1. So, at fixed A two roots merge at a point 
r = rtp > 0. The turning point approaches r = 0 in the limit A —> 0.) For later (preferred) times, they return towards the universal 
horizon, approach it asymptotically for v —> —oo, —>■ -|-oo, and are highly blueshifted. In addition, for r < M, there is a new v mode 
with negative norm for A > 0. It is indicated by a dashed green line in Fig. 4. It emerges from the singularity and approaches the 
universal horizon while closely following the positive Killing frequency characteristic after its turning point. (In fact this new v mode 
is directly related to the u modes emerging from the singularity in [2]: inside a universal horizon, u and v modes are swapped because 
of the vanishing of c at r = M.) Since some of the v modes originate from the singularity, and since the blueshift they experience is 
unbounded for r M~, the v part of the state will not obey Hadamard regularity conditions. This strongly indicates that the inner 
side of the universal horizon should be singular. This interesting question goes beyond the scope of the present paper. 
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FIG. 5. Characteristics crossing the infalling shell in the v — r,r plane. The Killing frequencies of the outgoing u modes and the 
incoming v modes is A = ±0.5A. The solid (dashed) lines represent characteristics for which the value of the Killing frequency 
A of the out-going u mode is positive (negative). The arrows indicate the future direction associated with the aether field. 
When tracing backwards the u-like characteristics associated with the Hawking quanta (A > 0) and their inside negative energy 
partners (A < 0), we see that they both originate from infalling u-like superluminal characteristics with a high and positive 
Killing frequency Ain. The v mode which emanates from the singularity (the dashed line) returns to it after having bounced at 
the center of the shell (not represented), closely following the characteristic of the v mode coming from infinity which hits the 
shell at the same value of r. 


4- The characteristics in the collapsing geometry 

In the internal region v < 4M, the characteristics are straight lines. Coming backwards in time from the outside 
region, the inside trajectories are fixed by the value of the inside frequency which is determined (as in the relativistic 
case), by continuity of the field tp across the mass shell; see Appendix A3 for details. As a result, the derivative i9r'0 
must be continuous across v = 4M. At the level of the characteristics (i.e., in the geometrical optic approximation), 
this implies that ky, the radial momentum at fixed v, is continuous along the shell. In terms of the inside and outside 
preferred momenta P'^{uj) and P“(A,r) evaluated at u = 4M“ and v = 4M+, respectively, the continuity condition 
gives 


■-M 


(A + P“(A,r)) =w + P“(a;). 


(27) 


This equation has two solutions, but only one is well behaved as the other one gives a trajectory along which the 
preferred time is not monotonic. A straightforward calculation using the dispersion relation Eq. (21) also shows that 
the sign of is preserved. It should be noted that Eq. (27) is the dispersive version of the relativistic equation 
|r/(r — 2M)|A = ui, which gives back Eq. (12) for r > 2M, cu ^ A, and when using u rather than r. 

It should be also emphasized that all outgoing u-like characteristics originate from inside the shell, as in the 
relativistic case. This is shown in Fig. 5. Therefore, thanks to the universal horizon, the state of the field inside the 
shell determines the state of the u modes. In this we avoid the problem discussed in [2], namely that in the absence 
of a universal horizon, the u modes of a superluminal field originate from the singularity at r = 0. As discussed in 
footnote 3, these modes still exist, but they are now trapped inside the universal horizon. 

Finally, we notice that the Killing frequency Ain of the incoming v modes which generate the outgoing u modes 
exiting the shell at Tc ~ M is very large. More precisely, when dealing with u characteristics with positive 12 
(i.e., modes with positive norm), irrespective of the sign of their Killing frequency A, the Killing frequency Ain is 
positive. A straightforward calculation (based on the continuity of ky applied to the v modes) shows that it scales as 
Ain « 3AM/ {ry - M). 

For completeness, we have also represented in Fig. 5 a couple of infalling v characteristics which enter the shell for 
0 < r < M. One comes from r = 0 (the dashed line), and one from r = oo (the solid line). They both reach the 
singularity after having bounced at r = 0 inside the shell. These characteristics, although interesting, play no role in 
the Hawking process. 
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C. Behavior of the WKB modes near the universal horizon 


To be able to compute the late time behavior of the Bolgoliubov coefficients, we need to further understand the 
properties of the stationary modes in the immediate vicinity of the universal horizon at r = M. For r > M, the two 
roots and remain finite as r —)■ M. As can be seen in Fig. 4, the associated trajectories smoothly cross the 

horizon. They thus play no role in the large uj limit. In fact they describe out modes. 

The two other roots and both diverge as r ^ M. Importantly, they both satisfy 


pin _ 


KM r f M\ 

± + ^ 1-> 


r — M M 


r ) 


(28) 


where the + sign applies to and the - sign to P^’^ ■ We have added a superscript in to emphasize that this 
behavior is relevant at early time t, just after having crossed the shell. The simple relation between P^ and P^’^ 
implies that for r —>■ M, the two WKB modes ifj\ and are also related to each other by flipping the sign of A. 
In the forthcoming discussion, to implement these points, we shall replace '4’'^'^ by and add a superscript in to 
the WKB modes ip^x- 

The appropriate character of this superscript can be understood as follows. Although the divergence in l/(r — M) 
in Eq. (28) resembles to what is found in the relativistic case, it has a very different nature due to the different 
relationship between r and the preferred coordinate X. This can be seen by looking at the validity of the WKB 
approximation for V’a ™ close to the universal horizon. Deviations from this approximation come from terms in 
{dxr)/irP), {dx{r - M))/{{r - M)P), and {dxP)/P^- Using 


T — M 

dx = - dr, (29) 

r 

we find that these three terms go to zero as r —> M. Therefore, close to the universal horizon, the WKB approximation 
of Eq. (26) becomes exact for ■ In iact, these modes behave as the dispersive in modes near a Killing horizon [21, 
22]. Namely, they have a positive norm for all values of A and, moreover, contain only positive values of P“. We 
recall that this is the key property which also characterizes the so-called Unruh modes [26, 31] for a relativistic field. 

These are strong indications that no stationary emission should occur close to the universal horizon, as the pair 
production mechanism rests on deviations from the WKB approximation. This is confirmed in the next subsection. 


D. Bogoliubov coefficients from the scattering on the shell 


We are now in a position to determine the scattering coefficients which govern the propagation across the null 
shell. Inside the shell, one has the in mode Along the shell, for v = AM~, it is a plane wave which behaves as 

^ exp [i{uj + P“(w))r]. After having crossed the shell, for r/M — 1 ^ 1, it may be expanded in terms of the 
four WKB modes (which form a complete basis) 


J.U, in 


dX 




w,in 

A 


Ouj,X 


/ IUAn\^ 

(V'-A ) 


+ ■Auj,\'4^x T (V'a’ ) 


(30) 


We are interested in the coefficients 'yui,\ and Sui,x which multiply the two modes with divergent wave vectors and 
opposite norms. It should be pointed out that the integral over A runs from —oo to oo. The other two coefficients 
A(j,a and B,^ x multiply the two modes which remain regular across the universal horizon in the {v,r) coordinates. 
They vanish in the limit w —)■ oo. 

The calculation of jui.x and is straightforward in the (v,r) coordinates; see Appendixes A 4 and A 5. For 
jA] < A, we find that their ratio decays as 


5uj,x 

lin,X 



X exp(-2MP“(w)), 


(31) 


where P'^(w) ^ xfujK in the present high frequency regime. Equation Eq. (31) is the main result of the present work. 
Its meaning is clear: at late time, corresponding to the emission close to the universal horizon and thus to very large 
values of P^ ^ A/{r/M — 1); see Eq. (28), the propagation across the shell induces no mode mixing between the 
inside in-mode </>“’“ and the high momentum WKB mode with negative norm ('f/’”'™)*’ irrespective of the value (and 
the sign) of A. As a result, outside the shell, the state of the field is stationary, and the vacuum with respect to the 
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annihilation operators associated with V’a’*" ^^r A G (— 00 , 00 ). ^ It thus correspond to the in vacuum as described 
in [ 21 , 22 ], 5 


E. Genericness of Eq. (31) 


In this subsection, we distance ourselves from the model we considered to see how the above results may be affected. 
We first consider a modihcation of the mass shell trajectory close to the universal horizon. From the calculation of 
Appendixes A 4 and A 5, the factor exp 2 Ma/Aw^ in Eq. (31) comes from the fact that the phase of the mode 

inside the mass shell is 6 *int ~ — wT « a; (r — u), while that of the mode outside the mass shell is 0out ~ KM/x, where 
X = (r/M) — 1. At hxed v, we find that the stationary phase condition applied to 0in ± Oout (the upper sign applies 
to 7 while the lower sign applies to S) gives back the large frequency limit of Eq. (27) with x real for 7 , while x is 
purely imaginary for 5, with a modulus Let us now consider an arbitrary shell trajectory close to the universal 

horizon. We define an affine parameter y along this trajectory. The possible saddle points are located where 




AM\ 

X ) 


= 0 , 


i.e.. 


So, the location of the saddle is 


dT AM dx 
± —= 0 - 


dy 


dy 


^ , „_^dx A 


(32) 


(33) 


(34) 


We get the same result as before, up to the factor —M^. Therefore, the ratio |5 cj.a/7i.j,a| is still suppressed by an 
exponential factor in My/Aw, with a coefficient depending on the velocity of the mass shell when it crosses r = M. 

We now consider a generalization of the dispersion relation Eq. (21) with higher-order terms. Specifically, we 
consider the dispersion relation 


~/A A 20-1) ■ 

i=o 


Close to the universal horizon, the divergent wave vectors follow 

P ~ ±Ajv X 


-1 

±Av XN-l 


(35) 


(36) 


As before, the coefficient 7 corresponds to ± = -I- in Eq. (36). The value of the saddle point is then real, and the 
exponential factor appearing in 7^7 has a unit modulus. Instead, for the coefficient S, corresponding to the minus 
sign in Eq. (36), the solutions of the saddle point equation are 


X 


* 


i^l+2L 

(a) ^ ” 


I G Z. 


(37) 


Taking only the saddle points with negative imaginary parts, we find that is suppressed by a factor which is 
exponentially large in Interestingly, when using the inside spatial wave number P“(a;) rather than the inside 

frequency w, the norm of the coefficient 5uj,x always decreases as exp (—MAP'^{w)) with A > 0, which means that it 
is the diverging character of P“(a;) which guarantees that its sign does not flip when crossing the shell. 

Similarly, the exponential factor suppressing (5^^ is mildly affected by a change in the metric and/or the form of 
the mther field, provided the inside wave vector remains smooth, whereas the outside one diverges as a power law for 
r —> ruH^ where ruH is the radius of the universal horizon. This should remain valid as long as there is no divergence 
(or cancellation) preventing us from defining preferred coordinates in which the dispersion relation takes the form 
of Eq. (35) close to the universal horizon. Indeed, the construction of Appendix A 1 can be easily extended to a 
generic space-time with a Killing vector 7 , endowed with a generic timelike, normalized aether field . 


^ This conclusion differs from that reported in [18]. We do not understand the procedure adopted there, which apparently implies that 
the leading term in Eq. (28) does not contribute to the ratio of Eq. (31), thereby giving rise to a steady thermal radiation governed by 
the surface gravity of the universal horizon. Instead, the saddle point evaluation of 5,^ y performed in Appendix A 5 establishes that the 
leading term of Eq. (28) gives the exponential damping in of Eq. (31). 

® To be complete, one should propagate backwards in time the inside field configurations, and verify that they correspond to vacuum 
u-like configurations for r —> 00 , t —>■ — 00 . To verify this, we computed the scattering coefficients encoding a change of the norm of the v 
modes when crossing the shell. We found that they also decrease exponentially in y/ujA for cc —>■ 00 . We also recall here that the Killing 
frequency of the v modes engendering a stationary u mode diverges as Aj^ ^ 3AM/ (vc — M), where Tc is the radius when the u mode 
exits the shell. 
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IV. CONCLUSIONS 

We computed the late time properties of the Hawking radiation in a Lorentz violating model of a black hole with 
a universal horizon. To identify the appropriate boundary conditions for the stationary modes of our dispersive field, 
we worked in the geometry describing a regular collapse, and assumed that the inside state of the field is vacuum at 
(ultra) high inside frequencies oj A. We then computed the overlap along the thin shell of the outwards propagating 
inside positive norm modes, and the outside stationary modes. In the limit where the shell is close to the universal 
horizon, we show that the overlap between modes of opposite norms decreases exponentially in the radial momentum 
P. This result comes from the peculiar behavior of the momentum when approaching the universal horizon with a 
fixed Killing frequency; see Eq. (28). Although this behavior was found in a specific model, we then argued that it 
will be found for generic (spherically symmetric) regular collapses and superluminal dispersion relations. 

As a result, irrespective of the model, at late time, the state of the outgoing field configurations is accurately 
described, for both positive and negative Killing frequencies, by the WKB modes with large positive momenta P 
(and a positive norm). In this we recover the standard characterization of outgoing configurations in their vacuum 
state in the near horizon geometry. Indeed, the condition to contain only positive momenta P prevails for both 
relativistic and dispersive fields in the vicinity of the Killing horizon. The present work, therefore, shows that this 
simple characterization still applies in the presence of a universal horizon. 

Once this is accepted, the calculation of the asymptotic flux is also standard, and shows that for large black holes the 
thermality and the stationarity of the Hawking radiation are, to a good approximation, both recovered. This suggests 
that the laws of black hole thermodynamics should also be robust against introducing high frequency dispersion. 

As a corollary of the divergence of the radial momentum on both sides of the universal horizon, noticing that 
the inside configurations are blueshifted (towards the future), and that they have no common past with the outside 
configurations, it seems that the field state cannot satisfy any regularity condition across the universal horizon. It 
would be interesting to study the space of the field states, and determine whether some dispersive extension of the 
Hadamard condition can be imposed on the universal horizon. In the negative case, it seems that the universal horizon 
will be replaced by a spacelike singularity. 
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Appendix A: Wave equation and Bogoliubov coefficients 

In this appendix, we give the general formulas and main steps in the derivation of the results presented in Section HI. 


1. Preferred coordinates 

The preferred coordinates {t,X) are defined by the followng four conditions 

• = ±dx at fixed f; 

• dy = ±dt at fixed X; 

• dyT < 0 along the shell trajectory; 

• dyX > 0 along the shell trajectory. 

These four conditions uniquely define t and X as 

{ V — r, V < 4M, 

v — r’^, v>AMAr>M, (Al) 

— {v — r ^), V > AM Ar < M, 
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and 


X = 


r, V < AM, 

r’lj, V > AM Ar > M, 
—r^j, V > AM Ar < M. 


In these expressions, r^, = r + M \n\j^ — 1\ is the tortoise coordinate built around the universal horizon. 


(A2) 


2. Wave equation and scalar product 

The action Eq. (17) has a 17(1) invariance under $ -A from which we derive the conserved current density 

(V,$) + C.C., (A3) 

where “c.c.” stands for the complex conjugate, satisfying 

= 0. (A4) 

As the wave equation Eq. (18) is linear, one easily shows that defines a conserved (indefinite) inner product in 
the following way. Considering two solutions $i and $2 of Eq. (18), we first define ($i,$ 2 ) by replacing by 
<i)J and $ by ^2 in Eq- (A3). The inner product of these two solutions is then defined by 

{^,,^2)r = J d^xnpJf^{A>,,A>2), (A5) 

where rip is the unit vector perpendicular to the 3-surface defined by t = cst, and r is a time coordinate. When 
considering the 3-surfaces defined by t = cst., the above overlap simplifies and gives the standard (Hamiltonian) 
conserved scalar product of Eq. (20). 


3. Matching conditions on the mass shell 

In order to compute the overlap of two modes defined on either side of the mass shell, we need the matching 
conditions to propagate the modes from the internal region to the external one and vice versa. As we now show, 
they appear naturally when considering the behavior of 77" = J^dpV across the shell. To see this, we first rewrite 
J" ($i,$ 2 ) as 


77 " ($ 1 , $ 2 ) =-i (^V" -f A>1 - (h"'^Vp$ 2 ) (Vph^"V,$t)^ 

- (4>* o $ 2 ). (A6) 

Inspecting Eq. (18) and requiring that the second term has no singularity which cannot be canceled by the first one, 
we find that the quantities 

• «■, 

• Vph'’‘^Vcr<i), and 

• y=5(V° + ^hO"V,Vph^"V,) 

are continuous across v = AM. Since the complex conjugate of a solution of Eq. (18) is still a solution, this applies 
to $ = as well as $ = $ 2 . Therefore, in evaluating Eq. (A6) one can evaluate 4>* and the operators acting on it 
on one side of the shell, v = AM — e, e —)■ 0, while <I >2 and the operators acting on it are evaluated on the other side 
V = AM + e. 
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4. Calculation of 'yui,\ 

Let us consider two radial modes known on different sides of the mass shell: is known for v < AM and $2 

for V > AM. The complete expression of the scalar product in the v,r coordinates is somewhat cumbersome, but it 
greatly simplifies in the relevant limit where 

• ipi has a large frequency |a;| ^ A; 

• ip 2 has a large wave vector \ky^ 2 \ ^ A, A. 

We have introduced the wave vector = drS at a fixed v. For the modes we are interested in ky ^2 = (A 2 + P 2 ) /x 
and uj are of the order A/a;^, where x = (r/M) — 1. Keeping only the leading terms in the inner product then gives 

(^’l, ^’2)„ ~ ^ y + dr)^tjj* + y “ {dv + drf V'l 


M 


+■4)1 1 - d'^ll)2 - {dy + dr) Ipl [ -1 dr422, 


M 


(A7) 


with relative corrections of order x. When choosing for ipi the in mode of frequency uj, and for '02 the stationary 
WKB mode of Eq. (26) with the large momentum given by Eq. (28), we get 


(^ 1 ,^ 2 )^ ~ 47rM f dx( 
Jx>0 \ 


p3 p2 \ p 

UJ I ^ UJ I I UJ \ 1^1 

A 2 ^A^ ^ I 


p3 p2 \ p 

dx[^±-^±^ + ^ 

Ax x"^ x^ 


ji^^4:iM{uj-X)^-iM(uj-\-P^) r 

exp fi -h (2A ± A) M In |a::| — M (uj + F0) x 


(A8) 


In this equation, as well as in the remainder of this appendix, the sign ± discriminates between 7 and 6; see below. In 
the large frequency limit, we evaluate this integral through a saddle point approximation. The possible saddle points 
are the values of x where 


_d 

dx 


T - M (uj + Puj)x] Ki — [ =F - Mujx = 0, 

\ X J dx \ X J 


i.e.. 


A 

±-. 

UJ 


(A9) 


(AIO) 


This is very similar to the saddle point condition applied to the Bogoliubov coefficients describing the scattering of 
plane waves on a uniformly accelerated mirror [32, 33]. 

The coefficient 7 ^^^ is defined for ± 0 ; > 0. Since the integral runs over a: > 0, we must choose the saddle point x* 
at 

^ (All) 


We get 


7 w.A ~ 


I ^iM 

27r \uj 


exp 


iM [iuj - AX - Pyj^ 2a/A jwj + - (2A ± A) In ( — ) =F A 


It is easily shown that, under these approximations, the following unitarity relation is satisfied: 

poo 

/ duj-il^^-iy,^y ~ 5(X-X'). 

JO 


(A12) 


(A13) 


This implies that the coefficients are suppressed in the limit w —> 00 . 
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5. Calculation of 

The calculation of 5ui^\ follows the same steps. The saddle point equation now is 

A 


- 

— 




(A14) 


To be able to deform the integration contour to include the saddle point, we must choose the solution in the half-plane 
where the exponential decreases, i.e., 


x*s = -sgn{uj)i. 


(A15) 


The exponential factor in the integral then gives a suppression factor 


exp ^2 a/A |w| -I- TT ^-sgn (w) A -I- -I- A^ ^ . 


(A16) 


In addition, to the order to which the calculation was performed, the prefactor vanishes. As the first relative 
corrections from neglected terms are of order 0{x*) = k/ |a;|), we get 


= O 


f y/Mk\ 


exp - 


[-2M^/k 


w 


(A17) 


Appendix B: Acceleration of the aether field 


The acceleration of the aether field is 


Using Eq. (13), this gives for v ^ AM 




7^7m = 


M{i 


(Bl) 

(B2) 


For completeness, we now show that in 1 -|- 1 dimensions a stationary universal horizon requires that the aether field 
has a nonvanishing acceleration, thereby generalizing what was found in de Sitter in [7]. We consider a stationary 
space-time with Killing vector AT^, endowed with a timelike ®ther field u^. The universal horizon is defined as the 
locus where = 0. (Notice that the Killing field must thus be spacelike on the universal horizon.) In particular, 

cannot be aligned with . Using the Killing equation, the variation of along the flow of is 

(K"u,) = (B3) 


If is freely falling, 7 ^ = 0 and is tangent to the hypersurfaces of constant In particular, it is tangent 

to the universal horizon. In 1 -|- I dimensions, since AT'' and cannot be aligned, is not a tangent vector to 
the universal horizon, which is thus not stationary. Models with a stationary universal horizon are thus in a different 
class than those studied in [ 22 ]. 

To see the combined effects of the dispersion and acceleration, we show in Fig. 6 the local value of the wave vector 
in the v,r coordinates, for the outgoing u mode, as a function of r. We compare three models with the same 
parameters, and for A = 10“^A. The blue, solid curve shows the result for the model of Section III. The green, dotted 
curve shows the relativistic case. The red, dashed one shows the result for a dispersive model with a nonaccelerated 
preferred frame chosen to coincide with the aether frame of Section III at r = 2M.® We see in Fig. 6 that the three 
models give very similar results for r < 2M. Close to r = 2M, the relativistic wave vector diverges, while the 
nonaccelerated dispersive model still closely follows the accelerated one. When r is further decreased, the predictions 
of the two models separate: the nonaccelerated one gives a finite wave vector at r = M while the accelerated one 
gives fc oc (r — 


It must be noted that this model is not well defined for r oo. The reason is that at r = 2M, we have u ■ dv = 1/2. A nonaccelerated 
vector field w which coincides with u at r = 2M must thus satisfy the two conditions tu • Ot = 1/2 and w • w = 1 at r = 2M. From the 
free-fall condition, these two properties extend in the whole domain where the preferred frame is defined. Since they are incompatible in 
Minkowski space, we deduce that the domain in which the preferred frame can be defined does not extend to r ^ oo. A straightforward 
calculation shows that it extends up to r = 8M/3. However, as this model is well defined close to and inside the Killing horizon, it can 
be used to see the qualitative differences between the nonaccelerated and accelerated cases. 
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FIG. 6. Comparison of three wave vectors k{r, A)/A as a function of a; = r/M — 1 (both in logarithmic scales) in the relativistic 
case (green, dotted line), for freely falling preferred frame (red, dashed line), and in the model of Section III (blue, solid line). 
The Killing frequency is A = 10“^A, and the Killing horizon is located at x = 1. One clearly sees the unbounded growth of 
the relativistic wave number. More importantly, one also sees that the two dispersive wave vectors behave in the same manner 
across the Killing horizon. Hence the acceleration of has a significant effect on k only when approaching the universal 
horizon. 


Appendix C: Hawking radiation in the presence of a universal horizon 


In Section III, we showed that the late time emission from the universal horizon is governed by Bogoliubov coefficients 
which are exponentially suppressed when the inside frequency w 3> A. This result was obtained using the WKB 
approximation of the stationary modes just outside the universal horizon. This approximation is trustworthy as we 
verified that the deviations from the WKB treatment go to zero when approaching the universal horizon. This implies 
that at late time the inside vacuum is adiabatically transferred across the shell. Hence, the u part of the field state can 
be accurately described by the WKB high (preferred) momentum mode V'a ™ both signs of A. In this we recover 
the situation described in [21, 22, 29]. Therefore, the nonadiabaticity that will be responsible for the asymptotic 
radiation will be found in the propagation from the universal horizon to spatial infinity. The value of the Bogoliubov 
coefficients should essentially come from the stationary scattering near the Killing horizon. Hence, we expect to get 
a nearly thermal spectrum governed by the surface gravity of the Killing horizon, and with deviations in agreement 
with those numerically computed in [34, 35]. 

To verify this conjecture, we numerically propagate the outgoing mode from a large value of r/2M down inside 
the trapped region to r —?> M+. This mode can be written in the limits r ^ M and r —oo as 

P'Pir) ~ rx+A.r^, 

•}—yoo 

Plir) (ci) 

9 —yM 

where the WKB modes are as described in Section HI. The coefficient A\ governs the grey body factor. In our 
(l+l)-dimensional model, we have verified that it plays no significant role. (We found that \A\\^ is bounded by 0.16.) 
Hence, as usual, the Hawking effect is essentially encoded in the mode mixing of u modes of opposite norms. 

To efficiently perform the numerical analysis, we regularized the metric and aether field. In practice we worked with 
a metric of the form 

ds^ = (1 — 2/(r)) dv"^ — 2dvdr, (C2) 

and a unit norm aether field 


= dy- f{r)dr. (C3) 

These expressions generalize the model of Section HI which is recovered for /(r) = M/r. The Killing horizon 
corresponds to /(r) = 1/2, and the universal horizon to f(r) = 1. We can then define the preferred coordinate X 
along the lines of Appendix A1. For the numerical integration of Eq. (19), it is appropriate to work with / expressed 
as a known function of A. A convenient choice is 

(l _ [A 



fir{X)) = i 


(C4) 
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FIG. 7. (Left panel) Characteristics in the v — x,xIxk coordinates, where x = r — ru and xk gives the location of the Killing 
horizon. Only the region outside the mass shell and the universal horizon is represented. The parameters are A/k ~ 1.1, 
Xo — 0.5, and A/A = 0.01. (Right panel) Amplitude of the nonadiabatic corrections in the present model (blue, solid line) 
and the one from Section III (orange, dashed line) for A/A = 0.01, as a function of the preferred coordinate X — Xk, where 
Xk denotes the position of the Killing horizon. In this example, the surface gravity is k ~ A. We thus verify that the norm 
of the coefficient /3 a is of the order of the maximal value of the nonadiabatic parameter ~ 10, as expected from the analysis of 
nonadiabaticity [24]. 


77 is a positive parameter which must be equal to 1 to have an asymptotically flat space at r —>■ c» and a universal 
horizon at r —>■ M, i.e., X —>■ —00. In our numerical simulations, we worked with 77 < I to avoid large numerical errors 
due to the divergence of the dispersive roots P“a close to the universal horizon. We then checked that the scattering 
coefficients become independent on 77 in the limit 77 —>■ I, as should be the case since the WKB approximations 
become exact on both sides. The advantage of this model is that the metric coefficients and the aether held converge 
exponentially to their asymptotic values so that the asymptotic modes become exact solutions for r —)■ cx), provided 
the decay rate of the exponentially decaying mode is small enough. ^ The characteristics outside the universal horizon 
r = ru are shown in Fig. 7 . They exhibit the main properties of Fig. 4 . In the right panel we show the parameter 
governing nonadiabatic corrections, {dxP\) IP\, as a function of the preferred coordinate X in the present model 
and the one of Section III. In the present model, they go to zero exponentially both for X ^ 00 and X ^ —00. In 
the model of Section III, they decay exponentially at X —>■ —00 but only polynomially at X -^00. The two models 
become equivalent, in the sense that the value of Px{f) follows the same law close to the Killing horizon when working 
with the same surface gravity. 

The held equation was integrated numerically using [ 36 ], and the same techniques as in [ 34 , 37 ]. The results are 
shown in Fig. 8. We obtain two important results. First, at hxed A and k, the effective temperature dehned by 




1 

eVn/r/L’ 


(C 5 ) 


becomes independent of the regulator 77 as 77 —>■ 1 “. Second at low frequencies A <C A, we get a Planckian spectrum, i.e., 
Tx = constant, with deviations from the Hawking temperature compatible with the results of [8, 34 ]. This establishes 
that the propagation between the two horizons does not alter the thermal character of the outgoing spectrum. 

To conclude this numerical analysis, we numerically verify that the WKB approximation becomes exact when 
approaching the universal horizon. To this end, we show in Fig. 9 the logarithm of the relative deviation between 
the numerical solution and the weighted sum of the WKB waves ax'tpx’''' + ■ As X —>• —00, we clearly see 

that the numerical values of the deviations decay following the rough estimation of the WKB corrections given by 


{dxPx) / {P\)^ ■ (The relatively important spread and the plateau for X < — 1.8 seem to be due to numerical errors. 

Indeed, as the wave vector becomes very large, typically of order 10 ^, even a relatively small error in its value gives 
important and rapidly oscillating errors.) 


Notice that an exponential convergence for X —oo is required to have a universal horizon at a finite value of r = rjj, such that 
df/dr{r = ru) 7 ^ 0 . 


7 
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FIG. 8. Plot of the effective temperature T\ of Eq. (C5) divided by the Hawking temperature as a function of A/k, where 
K is the surface gravity. The values of A/k are 1/2 (solid line) and 3/2 (dashed line). For the smallest value of A, we show 
the results for rj = 0.9 (blue line), 0.94 (cyan line), and 0.98 (magenta line). For the largest value of A, these three curves 
are undistinguishable up to numerical errors. This indicates that the limit of the regulator t; —>■ 1“ is well defined, which was 
checked using a larger range of values for rj £ (0.8,0.99). Moreover, when increasing the dispersive scale A, we see that Tx 
closely agrees with the Hawking value K/2n for a larger domain of Killing frequencies. 



FIG. 9. Deviations from the WKB approximation. The solid line shows the natural logarithm of | {dxP\) / (Px)^ \ as a function 
of the preferred coordinate X in units of 1/A, and the points show the logarithm of the relative difference between the solution 
computed numerically and the corresponding sum of the WKB modes, for A = 1, Xq = 0.5, rj = 0.99, and A = 0.1. The 
important spread seems to be due to numerical errors a result of the increase of the momentum as X —^ —oo. The Killing 
horizon is located at X = 0. 
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